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ABSTRACT 

We investigate the constraints on neutron star mass and radius in GS 1826-24 from models of lightcurves 
and spectral evolution of type I X-ray bursts. This source shows remarkable agreement with theoretical calcu- 
lations of burst energies, recurrence times, and lightcurves. We first exploit this agreement to set the overall 
luminosity scale of the observed bursts. When combined with a measured blackbody normalization, this leads 
to a distance and anisotropy independent measurement of the ratio between the redshift 1 + z and color correc- 
tion factor fc. We find \ +z= 1.19-1.28 for /c = 1.4-1.5. We then compare the evolution of the blackbody 
normalization with flux in the cooling tail of bursts with predictions from spectral models of Suleimanov et 
al. (201 lb). The observations are well described by the models at luminosities greater than about one third of 
the peak luminosity, with deviations emerging at luminosities below that. We show that this comparison leads 
to distance independent upper limits on R^o and neutron star mass of R^o 9.0-13.2 km andM < 1.2-1.7 Mq, 
respectively, for solar abundance of hydrogen at the photosphere and a range of metallicity and surface grav- 
ity. The radius limits are low in comparison to previous measurements. This may be indicative of a subsolar 
hydrogen fraction in the GS 1826-24 photosphere, or of larger color corrections than that predicted by spetral 

models. Our analysis also gives an upper limit on the distance to GS 1826-24 ofd < 4.0-5.5 kpc ' , where 
^i, is the degree of anisotropy of the burst emission. 

Subject headings: stars:neutron - X-rays:bursts - X-rays: individual (GS 1826-24) - X-rays: stars 



1. INTRODUCTION 

Thermonuclear flashes from accreting neutron stars, ob- 
served as type 1 X-ray bursts, in principle provide a way to 
constrain neutron star masses and radii (for a review see §4 of 
Lewin, van Paradijs & Taam 1993). The large observational 
catalogues of type I X-ray bursts now available (Galloway et 
al. 2008; hereafter G08) and new spectral models (Madej, 
Joss, & Rozanska 2004; Majczyna et al. 2005; Suleimanov, 
Poutanen, & Werner 2011b) have motivated fresh attempts 
to do this using photospheric radius expansion bursts (Ozel 
2006; Ozel, Guver, & Psaltis 2009; Guver et al. 2010a, 2010b; 
Steiner et al. 2010; Suleimanov et al. 2011a, 2011b; Ozel, 
Gould, & Giiver 2011). In this approach, the peak luminosity 
of the burst (specifically at the "touchdown" point when the 
photosphere returns to the neutron star surface) is related to 
the Eddington luminosity, and the normalization of the burst 
spectrum is related to the emitting area. If some informa- 
tion about the distance to the source is available, constraints 
on the neutron star mass and radius can be derived. These 
works have highlighted and spurred debate about some of the 
systematic errors that must be taken into account, such as un- 
certainty in identifying the moment at which the photosphere 
touches down (Galloway, Ozel, & Psaltis 2008; Steiner et 
al. 2010; Suleimanov et al. 2011a,b; Guver, Ozel & Psaltis 
2011), and differences in derived radii when using bursts at 
different accretion rates from the same source (Suleimanov et 
al. 2011a,b; Guver, Psaltis, & Ozel 2011). 

GS 1826-24 is a unique X-ray burster that shows remark- 
able agreement with theoretical models of recurrence times, 
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energetics, and lightcurves (Galloway et al. 2004; Heger et 
al. 2007; in 't Zand et al. 2009). The observed recurrence 
time (typically 3-5 hours) in a given epoch is the same from 
burst to burst to within a few minutes (Cocchi et al. 2001), 
and the burst lightcurves in a given epoch are very uniform 
(Galloway et al. 2004), implying the same conditions on the 
neutron star surface from burst to burst, and a regular limit 
cycle. Heger et al. (2007) compared the observed lightcurves 
to the theoretical models of Woosley et al. (2004). By choos- 
ing a model with approximately the same recurrence time as 
the data, and by varying the distance only (which scales the 
observed peak flux to match the model peak luminosity), the 
theoretical lightcurve fit most of the observed lightcurve well, 
except for deviations during the burst rise and at late times 
in the cooUng tail. This is of great interest because the long 
~ 100 second tails of these bursts are powered by the rp- 
process (Wallace & Woosley 1981), and offer a way to test the 
nuclear physics input, such as masses of proton rich heavy nu- 
clei and their reaction rates, both of which are usually highly 
uncertain (e.g. Schatz et al. 1998; Schatz 2006). Even more 
remarkably, in 't Zand et al. (2009) extracted the observed 
lightcurve out to more than 1000 seconds by combining mul- 
tiple bursts. The late time cooling observed matched the late 
time cooling in the theoretical model of Heger et al. (2007), 
which arises from heat initially conducted inwards to deeper 
layers that then emerges on long timescales. 

In this paper, we take the comparison between observations 
and theory for GS 1826-24 one step further. We first de- 
termine the constraints on neutron star mass and radius that 
can be derived from the lightcurve comparison carried out by 
Heger et al. (2007). The basic idea here is that even though 
the X-ray bursts from GS 1826-24 do not reach the Edding- 
ton Umit, we can still determine the intrinsic luminosity of 
the bursts by comparing with the lightcurve models. This re- 
places the touchdown measurement used for PRE bursts with 
a different condition which we use to constrain M and R for 
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GS 1826-24. In our analysis we take care to include the pos- 
sible anisotropy in burst and persistent emission and show 
how that could affect the mass and radius determination. We 
then compare the spectral evolution during the tail of the burst 
with spectral models. The good understanding of bursts from 
GS 1826-24 suggests that they could be a good testing ground 
for spectral models. We show that in the initial cooling phase 
following peak luminosity the spectral evolution agrees well 
with the models of Suleimanov et al. (2011b), and we de- 
rive the associated constraints on M and R. In both cases, 
we look for constraints that are independent of distance and 
emission anisotropics since neither are well-constrained for 
GS 1826-24. 

The outline of the paper is as follows. The data analysis is 
described in §2. In §3, we discuss the possible anisotropy of 
the burst and persistent emission and review calculations of 
the expected degree of anisotropy in the literature. In §4, we 
use the model lightcurve from Heger et al. (2007) to set the lu- 
minosity scale of the observed bursts and show that this gives 
a distance-independent relation between the redshift and color 
correction factor fc. Suleimanov et al. (2011a) argued that 
rather than using a single measurement of touchdown flux, the 
entire cooling track of the burst should be fit to spectral mod- 
els. We do this in §5, and show that even though the peak flux 
is below Eddington, the fits provide a constraint on the value 
of -Fkdd as well as the normalization of the spectrum. These 
two measurements translate into a distance independent up- 
per limit on Rao- We compare these two different constraints 
and discuss their implications in §6. 

2. DATA ANALYSIS 

We used data taken with the Proportional Counter Array 
(PC A; Jahoda et al. 1996) onboard the Rossi X-ray Timing 
Explorer (RXTE), from the catalogue of bursts detected over 
the mission lifetime (G08). Where not explicitly stated, the 
data analysis procedures are as in G08. Time-resolved spec- 
tra in the range 2-60 keV covering the burst duration were ex- 
tracted on intervals as short as 0.25 s during the burst rise and 
peak, with the bin size increasing step-wise into the burst tail 
to maintain roughly the same signal-to-noise level. A spec- 
trum taken from a 16-s interval prior to the burst was adopted 
as the background. 

We re-fit the spectra over the energy range 2.5-20 keV us- 
ing the revised PC A response matrices, vII.T] and adopted 
the recommended systematic error of 0.5%. The fitting was 
undertaken using XSPEC version 12. In order to accommo- 
date spectral bins with low count rates, we adopted Churazov 
weighting. We modelled the effects of interstellar absorption, 
using a multiplicative model component (wabs in XSpec), 
with the column density Nh frozen at 4 x 10^' cm"^ (e.g. 
in 't Zand et al. 1999). In the original analysis carried out 
by G08, the neutral absorption was determined separately for 
each burst, from the mean value obtained for spectral fits car- 
ried out with the Nh value free to vary. This has a negligible 
effect on the fluxes, but can introduce spurious burst-to-burst 
variations in the blackbody normalisation. 

The burst data used here has been corrected for "deadtime", 
a short period of inactivity in the detectors following the de- 
tection of a X-ray photon. There are however concerns re- 
garding the absolute flux calibration of the PCA associated 
with variations in the flux from the Crab nebula and the effec- 
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tive area of the PCA. We will show that such absolute uncer- 
tainties will not influence our derived constraints. 

3. ANISOTROPY OF THE BURST AND PERSISTENT 
EMISSION 

The possibility that the burst or persistent emission is not 
isotropic has been long discussed (e.g. Lapidus et al. 1985), 
but has not always been included in recent work using X-ray 
bursts to constrain neutron star mass and radius. For exam- 
ple, in Ozel (2006), Steiner et al. (2010), and Suleimanov 
et al. (2011a) the burst emission is assumed to be isotropic. 
Here, we review the expected size of the anisotropy. We fol- 
low Fujimoto (1988) and define an anisotropy parameter ^ by 
the relation A-nd^F^ = L between the observed flux F and the 
luminosity of the source L over the whole sky, where d is 
the distance to the source. When ^ < 1 (> 1), the radiation 
is beamed towards (away from) the observer We write the 
anisotropy factor for the burst and persistent emission as 
and respectively. 

Lapidus, Sunyaev & Titarchuk (1985) showed that if the 
accretion disk extends to the neutron star surface during the 
flash, it will intercept w 1/4 of the radiation from the burst, 
reflecting it preferentially along the disk axis. They provide 
the approximate expression 

Cb' = \ + \cosi\ (1) 

where / is the inclination angle (/ = 0° means the system is 
viewed face on, looking down the disk axis), which closely 
fits their more detailed results derived from solving the ra- 
diative transfer equations for a disk geometry (they found a 
maximum value of 1.39 rather than 1.5). The range of is 
from 0.5 (edge on) to 1.5 (face on), implying an uncertainty 
of a factor of 3 depending on inclination angle. 

The anisotropy factor for the persistent emission is perhaps 
even more uncertain than that for the burst flux, depending 
on the specific model for the inner accretion disk, boundary 
layer, and corona etc. Lapidus et al. (1985) and Fujimoto 
(1988) derive opposite behaviors for the factor as a func- 
tion of inclination. The model presented in Fujimoto (1988) 
for the persistent emission assumes that radiation from the 
boundary layer, which encircles the neutron star in a "belt" 
about its equator, is largely screened by the inflated inner part 
of the accretion disk and scattered preferentially in a direc- 
tion along the disk axis. In Lapidus et al. (1985), however, 
the inner part of the disk is assumed to be thin, and less than 
one half of the boundary layer radiation falls on the accretion 
disk and is re-scattered, again preferentially along the disk 
axis, while the remainder of the emission is beamed prefer- 
entially in the direction / = 90° (along the plane of the disk). 
This difference in their modelling of the inner accretion disk 
is made apparent by fact that, while Fujimoto (1988) predicts 
no radiation to be emitted in the / = 90° direction, Lapidus et 
al. (1985) find a substantial portion of the persistent emission 
will be beamed in that direction. The ratio £,p/^h varies by up 
to a factor of ^3 with inclination for both models, although 
while Fujimoto (1988) finds that the ratio is monotonically 
increasing with inclination, Lapidus et al. (1985) find the op- 
posite trend. We note that these two models do not consider 
the effects of general relativity on the trajectories of photons 
near the neutron star surface. However, Lapidus et al. (1985) 
show that when considering this effect, a substantially larger 
proportion (^ 28% for a ratio of the neutron star radius to the 
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Schwarzschild radius R/rs of 3) of radiation falls on the ac- 
cretion disk, further enhancing beaming along the disk axis. 

The definition of ^ is such that it always appears with dis- 
tance d in the combination S}l^d. Therefore, the uncertainty 
in anisotropy factor acts in the same way as an additional 
uncertainty in the distance to the source. For GS 1826-24, 
Homer et al. (1998) suggest a limit ; < 70° based on the 
low amplitude of the optical modulation at the orbital fre- 
quency, in which case equation (mi gives 0.84 < ^jj^' < 1.5, 

1/2 

or 0.85 < < 1.1. Therefore even if the distance to 
GS 1826-24 was perfectly known, anisotropy of the burst 
emission would represent an uncertainty of about ±15% in 
any quantity that depends on distance. For example, using 
spectral fits to determine R^o is subject to this uncertainty 
since the normalization of the spectrum depends on the solid 
angle R^/d^(_. Given these uncertainties, in this paper we 
look for constraints on M and R that are independent of dis- 
tance and anisotropy. 

4. COMPARISON BETWEEN THE OBSERVED AND 
MODEL BURST LIGHTCURVES 

We first ask what constraints on M and R arise from the 
comparison between the observed lightcurve and theoretical 
models of Heger et al. (2007). Photospheric radius expansion 
(PRE) bursts are often used in work to constrain neutron star 
properties from X-ray bursts, because the peak luminosity of 
the burst can then be taken to be the Eddington luminosity. 
This cannot be done for GS 1826-24 because the bursts do not 
show PRE, implying that they have a peak luminosity below 
Eddington. Instead, here we pursue the idea that the model 
lightcurves which fit the observed lightcurves so well tell us 
the peak luminosity of the bursts. 

Heger et al. (2007) selected from their models one that 
had a similar recurrence time to the observed bursts in 2000 
(the model had f,ecm = 3.9 hr as opposed to the observed 
frecm = 4.07 hr). They showed that, when the distance to 

1 /2 

GS 1826-24 (actually d) is chosen to make the predicted 
peak flux match the observed lightcurve, the theoretical and 
observed burst lightcurves show remarkable agreement. They 
considered fixed values of M and R, but the choice of those 
two parameters also changes the mapping between observed 
and model burst fluxes. Therefore, rather than vary distance 
alone, we find the value of the ratio of the observed flux to the 
model flux 

^=4-'f^)'-^ (2) 

that gives the best fit between the model and the 
data. Taking the peak values, /^modeipk = 1.29 x 10^^ erg 
cm"^ s"' (computed from the redshifted peak luminos- 
ity quoted in Heger et al. 2007, with R = 11.2km and 
z = 0.26) and Fobs.pk = 2.84 x 10"^ erg cm"^ s"', we find 



^ibsMnodel = 2.20 X 10-33. 

tion (|2|, we find 



Substituting this value into equa- 



^ = 10.9 km/6.0 kpc \^ 



model 



2.2 X 10- 



1/2 



(3) 



where we use the redshift assumed by Heger et al. (2007). 
Note that the model lightcurve is likely to be insensitive to the 
model gravity, so that the ratio /\,bs/fIiiodei does not depend 
sensitively on the M and R used in the model. For exam- 
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Fig. 1. — The average burst profiles with /recur = 4.07 hr compared with 
three separate fits of the mean theoretical model lightcurve from Heger et al. 
(2007) (model A3 which had a similar recun'ence time) is shown in the upper 
plot, with an inset showing only the first 30 seconds. The model has been fit 
to the data by varying the overall normalization, start time, and redshift. The 
solid (with red band) and dashed (with green band) lines represent the fits 
to the entire lightcurve and the first 30 seconds, respectively, and the dotted 
line (with blue band) represents the lightcurve fitted only by matching the 
peak fluxes with a fixed redshift of 1 + z = 1 .26. The bands show the range of 
luminosity variations from burst to burst in the theoretical model. The lower 
plot shows the difference between each of the fits and the average observed 
lightcurve, with a horizontal solid line at SF/F = for clarity. 



pie, the ignition column depth is weakly dependent on gravity 
in this burning regime (Bildsten 1998 derives yig„ cx g'^^"^). 
However, this is something that should be explored in further 
simulations. For now, we assume /VibsZ-finodei is a constant, 
and take equation (j3]l as a joint constraint on R and l+z- 

The theoretical uncertainty in Fmodei is at present un- 
known. The predicted lightcurves depend on the input nuclear 
physics, and prescription for convection and other mixing pro- 
cesses for example. These prescriptions vary from code to 
code, and currently only simulations from the Kepler code 
(Woosley et al. 2004) have been compared to the observations 
of GS 1826-24. Further simulations and comparisons are re- 
quired to determine what range of predicted peak fluxes still 
produce lightcurves with the correct shape to fit the data. For 
now, in order to put an error bar on the prefactor in equation 
we assume that the theoretical uncertainty in /^obs/^modei 
is ±10%, and keep in mind the fact that this number is uncer- 
tain. 

This raises the point that rather than use the peak flux 
only, we could also fit the entire lightcurve. In that case 
there is an extra parameter, the redshift l+z which stretches 
the lightcurve in time. In principle, this provides a con- 
straint on l + z. In practice, however, we find that the value 
of 1 + z obtained in the fit is sensitive to how much of the 
lightcurve is included in the fit. For example, fitting the entire 
lightcurve (until about 130s after the peak) we find best-fit 
values 1 + z = 1.44, Fobs/fmodet = 2.10 X 10-^^. If we fit the 
first 30 seconds only, which includes only the initial decline 
after the peak rather than the whole tail, we get a best fit of 
1 +z = 1.32 and Fobs/-Fmodef = 2.17 x lO-^^. We show in Fig- 
ure [T| the separate fits to the entire lightcurve and the first 30 
seconds, and we also include the model lightcurve fitted only 
by matching the peak fluxes, with the value for the redshift 
of 1 + z = 1.26, as assumed by Heger et al. (2007) . We see 
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that while the redshift is sensitive to the details of the fitting, 
the normalization Fobs/^modei is well-determined. Therefore, 
here we use the normalization, but leave fits to the shape of 
the entire lightcurve to future work when a greater number of 
simulations are available. 

A second constraint comes from spectral fitting. Fitting the 
observed burst spectrum with a blackbody gives the total flux 
Fao, color temperature oo, and the blackbody normalization 



K = 



Ri 



d'fc 



•4?i 



(4) 



The color correction factor = r^ /Teff takes into account the 
hardening of the burst spectrum compared to a blackbody at 
the same effective temperature T^ff. We plot A' as a function of 
time in Figurel2]for two average burst profiles, for recurrence 
times 5.74 ancM-.O? hr respectively. Following the burst rise, 
which lasts for approximately 5 seconds, the normalization 
levels off until w 60 seconds into the burst, when the nor- 
malization drops dramatically over 100 seconds to only about 
25% of its original value. We discuss the variation of K in 
the tail, and the difference in K for the two different recur- 
rence times in the next section. For our purposes here we 
find the mean value of K for the f,ecm = 4.07 hr profile dur- 
ing the period following the peak where it is constant, giving 
K= 1 10 ± 2 (km/ 10 kpc)^. If we take into account the dead- 
time correction near the burst peak 6%), the value we find 
is consistent with the more detailed analysis of blackbody nor- 
malization in these bursts carried out by Galloway & Lampe 
(2011). 

Dividing equations (j2]) and Q, R/d and drop out, giving 



1+Z 



= 1.17 



Fobs 



model 



1/4 



(5) 



K 



110(km/10kpc)2 



-1/4 



2.2 X 10- 



1/4 
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We show model calculations of fc from Suleimanov et 
al. (201 lb) in the next section which typically have fc ~ 1 .4— 
1.5 during the phase where K is relatively constant. For 
fc = 1.4,1.5 we get l-l-z= 1.19,1.28. Note that as well 
as being independent of d and the value of 1 + z de- 
termined in this way is not very sensitive to the values of 
K and f,nodei/^obs (proportional to the 1/4 power of each). 
For example, introducing an uncertainty in /^modei/^obs of 
±10% gives a prefactor in equation (j9]l of 1.17 ±0.03 or 
l+z= 1.28 ±0.03 (/c/1.5). 

There is one more constraint which comes from the agree- 
ment between the measured and model recurrence times, 
which effectively measures the local mass accretion rate m 
onto the star. We define m to be the rest mass accretion rate 
at the stellar surface. Then the accretion flux as observed at 
infinity is 



ch 



(7) 



z, 



(8) 



Dividing equations (|2]l and (j7]l gives the observed quantity 

Fx Fjnodel _ f 

mc^ Fobs \£.p 

a direct measure of redshift, independent of distance, but 
dependent on the anisotropy parameter ratio £,p/£,b- The 
accretion rate in model A3 of Heger et al. (2007) was 



."3.74 hr recurrence time bursts □ 
4.07 tir recurrence time bursts O 




Time (s) 

Fig. 2. — Blackbody normalization for average burst profiles with 
cm =5.74 (black squares) and 4.07 hr (red diamonds). 



m = 7980 g cm ^ s ' , and the measured persistent flux in 2000 
was Fx = 2.91 ± 0.03 x 10"^ erg cm"^ s"', giving 



z = 0.18±0.02 



(9) 



An alternative way to derive this result is to match the theoret- 
ical a value, the ratio of persistent fluence between bursts to 
burst fluence, to the observed value. The models of Heger 
et al. (2007) with recurrence time of 4 hours have a theo- 
retical value amodel = '^Mc^Zmodel/EbuiA^.b/^,,) = 55 (£.b/ip), 

where finuc is the burst energy, AM the ignition mass, and 
a redshift Zmodei = 0.26. The observed a at the same recur- 
rence time is aobs ~ 37 (Fig. 2 of Heger et al. 2007), giving 
z = Zmodei(aobs/amodei) = 0. 17(^^,/C/,), in agreement with the 
value in equation (j9]). 

If the anisotropy parameters were known, equations (19]), (j2]i 
and (pl uniquely determine the three quantities 1 + z, /c and 
However, as noted in §2, the anisotropy parameters are 
not well-constrained. In particular, equation ^ for the red- 
shift is not constraining once the large uncertamty in £,p/£,i, is 
taken into account. The main result of this section is therefore 
the relation between 1 + z and fc in equation (j5]l, since it is 
independent of d and £_h- 

In the next section, we constrain the Eddington flux by fit- 
ting the burst cooling tracks to spectral models. We can use 
the model lightcurve to say something about the expected 
value of /^Edd, the observed flux which corresponds to the 
Eddington flux at the surface of the star In the model the 
Eddington flux locally is F^dd = cg/n = 0.882 x IO^^^h erg 

Cm-2 S"' for X = 0.7, giving ^model.pk/^model.Edd = 1 •46/^14. 

Scaling the observed peak flux, the Eddington flux as ob- 
served at infinity should be 1.95 x 10-'^gi4(1.7/l -l-X) erg 
cm-2 s-i, or /^Edd/(10-'* erg cm'^ s"') = 1.95,3.89,7.79 for 
logiog= 14.0, 14.3, 14.6. 

5. COMPARISON WITH SPECTRAL MODELS 

We now turn to fitting theoretical calculations of the color- 
correction factor fc to the data. First we describe the fit- 
ting procedure and results (§5.1) and then the constraints on 
neutron star parameters, in particular an upper limit on R^o 
(§5.2). In §5.3 we discuss the variation of K with accretion 
rate (Fig. 2) in the context of the spectral models. 
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Fig . 3 . — Best fits to the theoretical fc—F/ FEdd curves for a range of compositions. The crosses represent the all the data points from the burst peak onwards, and 
those in red representing only the first 35 seconds after the peak. Two vertical dotted lines represent the fluxes at ? = 10, 35 s after the burst peak. The compositions 
are solar H/He abundance with solar metallicity (diamonds connected by purple lines) or 1% solar metallicity (triangles, blue lines), pure H (squares, red line) 
and pure He (x symbols, green line). Dotted, dashed and solid lines represent surface gravities of log[Q(g) = 14.0, 14.3 and 14.6, respectively. 



5.1. Fit for A and F^^^ 

Suleimanov et al. (2011b) calculated fc as a function of 
F I fedd for a range of surface gravities and atmospheric com- 
positions, and discussed how these models could be applied 
to data. We follow their analysis, and fit the theoretical f^- 
F I FEdd curves to the observed relation between /T"'/'* and flux 
F . The fitting parameters are Fgdd and A = K''^!"' / f^. Compar- 
ing with equation (|4]i, we see that 

If fc was a constant independent of flux, and K was constant 
in the cooling tail of the burst, fitting for A would be equiva- 
lent to using the measured normalization K and a value of fc 
to extract Roo/d. Instead, here we are using the entire cooling 
track to obtain A. In addition, by fitting the shape of the cool- 
ing track we can obtain the overall flux scale F-em even though 
the burst itself does not reach Eddington luminosity. 

We start by fitting the data from GS 1826-24 with recur- 
rence time 5.74 hr to the different models from Suleimanov 
et al. (201 lb). Below we use the fits to obtain an upper limit 
on R^, which motivates us to start with the bursts with the 
largest value of K and therefore larger R^o values. At the 
end of this section, we discuss whether it is possible to in- 
clude the 4.07 hr recurrence time bursts which have smaller 
values of K (Fig. 2) in a consistent picture. When fitting, for 
simplicity we calculate based on comparing /T"'/'* and fc, 
and do not include the errors in the flux measurement. This 
seems reasonable because the observational error in the over- 
all flux scale is ss 5F /\/N, where the individual flux error 
5F « 10"^ erg cm~^ s"', smaller than the overall uncertainty 
in the parameter fedd that we obtain from our fits. 

Suleimanov et al. (2011b) calculate spectral models for 
pure H and pure He atmospheres, and solar H/He fractions 
with different metallicity. Based on the lightcurve compar- 
ison and energetics. Galloway et al. (2004) and Heger et 
al. (2007) concluded that the accreted layer has solar metal- 
licity and a substantial amount of hydrogen (a solar H/He ra- 
tio in their models). Here we fit to the full range of models 
from Suleimanov et al. (2011b) to investigate how changing 
composition affects our derived limits on neutron star param- 



eters. Also, the photospheric abundances could be different 
from the abundances near the base of the layer where the X- 
ray burst ignites. For example, the metallicity in the burn- 
ing layer could be enhanced by partially burned fuel left over 
from a previous burst. The hydrogen fraction in the accreted 
material could be lower than solar, and so it is useful to con- 
sider the pure He limit as a limiting case when the hydrogen 
fraction at the photosphere is reduced. 

Figure [3] shows example fits to the 5.74 hr recurrence time 
bursts. By varying A and Fedd, we are able to obtain good 
fits for fluxes down to approximately 1/3 of the peak flux. 
At lower fluxes, the behavior of the model and observations 
is qualitatively similar, in that fc rises rapidly at low fluxes, 
but the detailed behavior does not match the models. At late 
times or low fluxes, K~^l'^ rises more rapidly than predicted. 
We therefore confine the fit to the initial part of the cooling 
tail and use it to derive A and /^Edd- To do this, we fit Afpt sec- 
onds of data starting at the time of peak flux. The time Affit is 
chosen so that as much data is included in the fit as possible 
while still giving an adequate fit (with the late time data ex- 
cluded, we find reduced values in the range 0.23-0.47 for 
the models listed in Table 1). For all except the solar metallic- 
ity models we take Affit = 35 s, corresponding to fluxes down 
to w 1/3 of the peak flux. The solar metallicity models be- 
gin to deviate from the data after Affit = 10 s (about 1/2 of the 
peak flux) because of the dip in fc at low fluxes. 

The results of the fits for different spectral models are 
listed in Table [T| We used Markov Chain Monte Carlo 
implemented with the Metropolis-Hastings algorithm (Gre- 
gory 2005, Chap. 12) to sample the parameter space and find 
the distributions for the fitting parameters, A and /^Edd- Those 
distributions were then each fitted by a Gaussian profile in 
order to derive their respective central values and Icr uncer- 
tainties. In certain cases, as noted in Table[T] the distributions 
were not well described by a single Gaussian profile, due to 
the presence of more than one peak. In those cases, we fit a 
Gaussian profile to the peak at the lowest values of F^a and 
A. Note that these parameters are correlated since an increase 
(decrease) in A, which moves the model curves upwards 
(downwards) with respect to the observations, can be offset 
by a corresponding increase (decrease) in /^Edd which moves 
the curves rightwards (leftwards). The range of /^Edd for all the 
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fitted models is from 4.1 to 7.4 x 10"^ erg cm"- s"', which lies 
within the range of pEdd from the Heger et al. (2007) models 
used in section §4. Excluding the 0. IZq, logjQg = 14.0 model, 
the range of F^dd is relatively narrow; 4.1 to 5.9 x 10"^ erg 
cm"^ s"'. 

The behavior at low fluxes is shown in more detail in the 
second panel of Figure |3] which has a logarithmic flux axis. 
The slope of the increase m fc with decreasing flux is steeper 
in the data than in the models for low metallicity models. For 
solar metallicity models, the slopes are similar This would 
enable a good fit of the whole data set to those models, par- 
ticularly at low fluxes, but only if Fedd is in the range 15- 
25 X 10"'^ erg cm"^ s"'. In fact, even with the restriction of 
using only 10 seconds of data following the burst peak we 
found other adequate fits, as separate local minima, in that 
range of /^Edd- This is much larger than expected and as can be 
seen from the relations derived below, would give very small 
limits on Roo, and so we do not consider these fits further 

5.2. Upper limit on distance and Roo 

A measurement of A and F-em translates into typically two 
values forM and R as follows. When the flux at the surface of 
the star is at the local Eddington flux eg/ k,, the observed flux 
is 



GMc 1 1 

nd^ 1+Z 



(11) 



where we take the opacity to be k = 0.2 cm^ g"' (1 +X), 
as used in Suleimanov et al. (201 lb). Following Steiner et 
al. (2010) we define the quantities 



nd , 



-FEAdA'i 



(l-«) 



R 



7^ 



'tA^FEdd (m/2)(1-m)3/2 



(12) 



(13) 



where u = 2GM/Rc^. These definitions differ slightly from 
those of Steiner et al. (2010) in that they include the 
anisotropy parameter ^h- Then 



«=i±l(l-8a)'/2 



R = a"f(l — u) 



1/2 



(14) 
(15) 



To calculate a and 7 from A and ^Edd obtained from the fits, 
we require distance d and composition X. A given a and 7 
then give two solutions for M and R. 

We treat the A and F^dd values as given independently of the 
derived M and R. In fact the derived A and F^dd values depend 
on the gravity assumed for the spectral models. The color 
correction fc decreases with increasing gravity in the models 
of Suleimanov et al. (2011b) (see their Fig. 5). This implies 
an increasing A with gravity, which we find in our results for 
the Z = O.OIZ0 models. A similar trend is not seen in the 
Z = O.IZ0, Zq models. The color correction still decreases 
with gravity, but the shapes of the models change in such a 
way as to favor smaller values of ^Edd, which, given the corre- 
lation between the two fitting parameters, leads to values for 
A that are smaller than expected for the higher gravities. The 
value of F-Edd also increases with gravity for the Z = O.OIZ0 
models, because the slope of fc W\\hF /F^dd steepens with in- 
creasing gravity, requiring a larger value of F^dd to agree with 
observed K~^l'^-F slope. 



Overall, we found that the fits are sensitive to the detailed 
shape of the atmosphere models. This could be due to some 
small irregularities in the model slopes attributable to the 
coarseness of the flux grid on which the color corrections 
are evaluated, and the limited range of fluxes spanned by the 
bursts we are analyzing. 

Table [T] shows that for a given metallicity, varying the sur- 
face gravity changes changing gravity from logjgg = 14.0 to 
14.6 changes A by up to « 3% for low metallicity, and ^Edd 
by up to w 20%. The resulting changes in the limit on R^o are 
sa 30% for low metallicity models. This gives a measure of 
the error introduced by not carrying out a self-consistent fit in 
which the gravity of the spectral model used and derived M 
and R are consistent. 

Equation ( [T2] i shows that real-valued solutions for u require 

< 1/8 (Steiner et al. 2010). As emphasized by Suleimanov 
et al. (201 la, 201 lb), this gives an upper limit on the distance 
for which solutions are possible. 



\c' 



1 



8 K FEddA2 



= 5.6 kpc 



/ -pEdd.- 



1.2/ V 4.0 



1.7 



(16) 



, (17) 



where As = A/10« and F^dd.-i = pEdd/ilQ-^ erg cm'^ s"'). It 
also provides a limit on Roo = R{1 +z). To see this, note that 
the neutron star radius is 



Roo=R{l+z) = aj 



= 12.0 km a 



1.2 



-4 



fkdd,-8 y' /i+^y 
4.0 / V 1.7 / 



(18) 
■(19) 



where the definition of 7 from equation ( 1 3 1 was used and 
(I-m)'/^ was substituted with (1 +z)"'. An upper limit on Roo 
is obtained by setting a = 1/8 in equation ( 19 1. 

The upper limits on £,1/^'^ and Roo are given in Table [1] To 



calculate them we use equations ( 16 1 and ( |T9] l with 95% Tower 
limits on the quantities A^/^Edd and A'^FEdd derived from our 
fits. A slightly different procedure is used for the cases where 
the fits yielded multiple minima. To derive the most con- 

1 /2 

servative upper limits on d and Roo, we consider only the 
X^ local minimum located at the lowest value of /^Edd and A, 
manifested as a distinct, Gaussian-like peak in the respective 
distributions for the quantities A^/^Edd and A'^/^Edd- Consider- 
ing only the part of the Gaussian-like distribution lying below 
the peak value, we find the 90% lower limits for A^/^Edd and 
A'^pEdd- This is equivalent to taking the 95% lower limit of the 
whole peak, but has the advantage of allowing us to isolate the 
X^ minimum of interest from the rest of the distribution. As 
a check, we applied this procedure to model fits showing a 
single x^ minimum, and found very small differences (< 1%) 
in the derived upper limits when compared to those found by 
considering the entire distributions. 

An upper limit on Roo implies an upper limit on the neu- 
tron star mass M^ax = c^Roo/i^^'^G (at that mass the radius is 
P^oo/%/3), also given in Table [l] Note that the upper limits 

on Roo and d are correlated. Since ^y^c/um = 7A^/8 (compare 
eqs. p3) and p6)), we can rewrite equation ([T9| as 



Roo < 12.0 km 




(20) 
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TABLE 1 

Fits of spectral models to the burst cooling tail 



Composition 






A 




xLL,ced(^-'^-/) 


(kpc) 




Upper limits 








(s) 


(lO**) 


(10"^ erg cm"2 s"') 




Roo (km) Mmax(A/0)'' 


Rm>Mq {kmf 


PureH 


14.3 


35 


1.144±0.002 


5.15±0.20 


0.27 (59) 


4.3 


10.2 


1.3 


8.2 


Pure He 


14.3 


35 


1.186±0.002 


4.14±0.09 


0.26 (59) 


9.7 


21.5 


2.8 


19.8 


O.OIZq'' 


14.0 


35 


1.138±0.002 


4.71±0.14 


0.24 (59) 


5.5 


13.2 


1.7 


11.3 


O.OIZq'' 


14.3 


35 


1.154±0.002 


5.08±0.15 


0.23 (59) 


5.0 


11.6 


1.5 


9.7 


O.OIZq'' 


14.6 


35 


1.167±0.002 


5.43±0.19 


0.26 (59) 


4.6 


10.5 


1.4 


8.4 


O.IZ©'' 


14.0 


35 


1.170±0.003 


7.35±0.66 


0.47 (59) 


4.0 


9.0 


1.2 


6.8 


O.IZq'' 


14.3 


35 


1.168±0.004 


5.80±0.41 


0.45 (59) 


4.5 


10.2 


1.3 


8.1 


O.lZ©'' 


14.6 


35 


1.174±0.003 


5.75±0.31 


0.36 (59) 


4.4 


9.9 


1.3 


7.8 


Zq<'.'= 


14.0 


10 


1.178±0.020 


5.92±1.22 


0.36 (31) 


5.3 


12.4 


1.6 


10.6 


7 d,e 


14.3 


10 


1.159±0.017 


4.81 ±1.08 


0.29 (31) 


5.6 


13.2 


1.7 


11.3 


7 d,e 
Zq 


14.6 


10 


1.164±0.011 


4.84±0.52 


0.25 (31) 


5.4 


12.6 


1.6 


10.7 



" We fit to data from the time of peak luminosity until Afflt seconds later 

The maximum neutron star mass consistent with the upper limit on Roc = R(l +z). 

The upper limit on radius assuming a lower limit on mass M > 1 M© . 
'' The composition is solar H/He abundance plus the indicated proportion of solar metallicity. 

" These fits yielded more than one local minima. Here, we report only the minima located at the lowest values of Fedd and A. See text for more details. 



a larger distance limit allows larger radii. 

For solar abundance of hydrogen at the photosphere, we 
find < 4.0-5.6 km and R^^ < 9.0-13.2 km. This rep- 
resents quite stringent limits on the neutron star mass and ra- 
dius. For this range of Roo, the maximum neutron star mass 
is in the range 1.2 to 1.7 Mq. If we consider a lower mass 
limit of 1 Mq, the neutron star radius must be smaller than 
R{1 Mq) = 6.8-11 .3 km (the individual values for each model 
are given in Table [TJ. 

5.3. Variation ofK with accretion rate 

Figure |2] shows that the bursts with recurrence times of 
4.07 hr have significantly smaller values of K than the 5.74 
hr bursts, by w 20% (see Galloway & Lampe 201 1 for a de- 
tailed discussion of the variation of K in the sample of bursts 
from GS 1826-24). Variations in K between bursts has been 
seen in other sources. For example, Damen et al. (1989) found 
that the blackbody temperature (evaluated at a fixed flux level) 
depended on burst duration. They suggested that variations 
in chemical composition at the photosphere and the resulting 
changes in color correction might explain the changing black- 
body temperature (and therefore normalization). 

We investigate two possible composition variations: chang- 
ing metallicity with solar H/He abundance, and changing the 
hydrogen fraction. First, we consider solar H/He abundance 
and changing metallicity. Suleimanov et al. (201 lb) show that 
fc drops with increasing metallicity. Therefore we fit the so- 
lar metallicity model to the 5.74 hr bursts to determine values 
of A and pEdd (as given in Table 1). These values are then 
used to compare a low metallicity model to the 4.07 hr data. 
This comparison is shown in the top panel of Figure H] The 
low metallicity model lies below the 4.07 hr data, snowing 
that the difference in K cannot be explained by a decrease in 
metallicity from solar to a fraction of solar. 

Second, we consider a change in hydrogen fraction at the 
photosphere. The lower panel of Figure |4] shows the pure He 
atmosphere fit for the 5.74 hour bursts (see Table 1), and a low 
metallicity solar H/He abundance model for the 4.07 hr bursts 
in which we use the same value of A determined by the pure 
He atmosphere fit, but decrease the derived /^Edd by a factor of 
1 +X = 1 .7 to account for the difference in Eddington flux with 
composition. This plot shows that the change in fc in going 




Fig. 4. — Top panel: The solar metallicity fit (log^Q g = 14.3, linked blue di- 
amonds) that reproduces the first part of the cooling track for the /^recur — 

5.74 

hi' bursts (blue crosses, and delimited from the data exluded from the fit by a 
dotted vertical blue line) is plotted together with the data for the ?iccui = 4.07 
hi' bursts (red crosses) and the low metallicity Z = O.OIZq model (linked red 
triangles) at the same A and fgdd th£ solar metallicity model. Bottom 
panel: A pure He spectral model (linked blue X symbols) fit to the 5.74 hr 
bursts (blue crosses) and a low metallicity solar H/He composition model 
(linked red triangles) at the same A and F^dd adjusted for the different hydro- 
gen fraction. In both panels, the respective symbols show the points where 
the atmospheric models were calculated. 
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Fig. 5. — Summary of distance-independent constraints in the neutron star 
mass-radius plane. The dashed curves are lines of constant surface gravity 
log[Q(g) = 14.0, 14.3, 14.6 (bottom to top), values at which the spectral mod- 
els were evaluated. In green, we show the redshift from eq. |5} for fc = 1.4- 
1.5 and an assumed 10% uncertainty in Fobs/^model. The squares (dark blue), 
diamonds (light blue) and triangles (purple) represent the upper limits on i?oo 
computed from fits to the solar H/He abundance models with 0.0 IZq , 0. \Zq 
and Zq metallicities, respectively, each at a specific surface gravity. The up- 
per limit on Roo for the pure Helium atmosphere model (log = 14.3) is also 
shown as a black asterisk. Two constant Roo curves are plotted as dotted lines 
for the highest and lowest values found within solar H/He abundance models. 
The region hashed in black represents what is allowed by the combination 
of the constraints derived from the fit to the burst lightcurve and spectral fits 
to solar H/He abundance models. These constraints are independent of the 
source distance and anisotropy parameters 5;,, 5,,. The region in red repre- 
sents the mass-radius relation derived by Steiner et al. (2010) (based on the 
'■ph 2> R assumption), with the Icr and 2a regions delimited by solid and 
dot-dashed lines, respectively. 

from pure He to solar H composition is enough to account for 
the variation in K observed. However, the solar composition 
model does not match the 4.07 hr data in terms of location 
on the F /Fem axis. Another way to say it is that if we fit the 
4.07 hr data with a solar composition model, the required /^Edd 
would be larger than for the 5.74 hr data, instead of being a 
factor 1 +X times smaller, as is required for simultaneous fits. 
Furthermore, we see in the lower panel of Figure |4] that re- 
ducing the derived Fedd by a factor of 1 +X = 1 .7 for the 4.07 
hr bursts implies that the peak flux for those bursts exceeds 
the Eddington limit, which is known not be the case. There- 
fore a consistent explanation of the variation in K in terms of 
changing H fraction at the photosphere is not possible. 

6. SUMMARY AND DISCUSSION 

We have compared lightcurve and spectral models with ob- 
servations of Type I X-ray bursts from GS 1826-24. Here 
we summarize the main conclusions and discuss our results 
further. 

A general point is that anisotropy in the burst emission en- 
ters as an additional uncertainty in any derived quantity that 
depends on distance. Since it changes the relation between 
the source luminosity and observed flux, the anisotropy pa- 
rameter (defined in §3) always enters in combination with 

1 /2 

distance as d. Even in cases where distance to a source 
can be accurately determined, the anisotropy introduces an 
effective uncertainty of up to a factor of 20-30%. Anisotropy 
could be a smaller effect for PRE bursts if the inner disk is 
disrupted during the burst and intercepts a smaller amount of 
radiation than a disk extending all the way to the stellar sur- 
face. Nonetheless, it remains a source of systematic error on 
derived neutron star radii that needs to be investigated further. 
For GS 1826-24, the Hmit / < 70° from Homer et al. (1998) 



-1/2 

gives ' = 0.9-1.2. Given this uncertainty and the fact that 
the distance to GS 1826-24 is not well constrained, we fo- 
cused on deriving limits on M and R that are independent of 
distance and anisotropy. 

The first of these constraints comes from using the model 
lightcurve from Heger et al. (2007) to fix the overall luminos- 
ity scale of the observed bursts. We showed that this leads 
to a distance and anisotropy independent relation between the 
redshift 1 -l-z and color correction factor fc (eq. ||5)) that de- 
pends weakly on the measured normalization K and the ratio 
of observed and model peak fluxes. For a color correction be- 
tween 1.4 and 1.5, which spans the range of values in Fig. 2 
of Suleimanov et al. (201 1) for example, the inferred redshift 
is between z = 0.19 and 0.28. 

The second constraint comes from comparing the spec- 
tral evolution during the cooling tail with the spectral mod- 
els of Suleimanov et al. (2011b), which determines the Ed- 
dington flux /^Edd and the quantity A = K~^l'^ / fc- As noted 
by Suleimanov et al. (2011b), for a given set of measured 
^^Edd^ A parameters, there is an upper limit to the distance of 
the source beyond which there is no solution for M and R. 
We point out here that measuring A and F-em also places an 
upper limit on R^o = R{l+z) (and therefore upper limits on 
M and R for a given source). This limit is independent of 
distance and anisotropy and depends only on the measured 
values of A and F^m and the surface hydrogen fraction. For 
GS 1826-24, atmospheric models with solar hydrogen frac- 
tions give Rao < 9.0-13.2 km (Table[T]l which implies a neu- 
tron star mass M < 1.2-1.7 Mq and R < 6.8-11.3 km assum- 
ing a lower mass limit of 1 Mq. The corresponding distance 

Hmits ai-e d < 4.0-5.6 kpc 

Uncertainties associated with absolute flux calibration do 
not affect our results; they are equivalent to an incorrect mea- 
surement of the distance to the source, which our constraints 
are independent of. 

The constraints on M and R are summarized in Figure l5] 
We show the upper limits on R^o from Table[T]for all the soim: 
hydrogen composition models each plotted at the respective 
surface gravity and the pure Helium model with log^ = 14.3, 
and the redshift range 1 -l-z = 1 .16-1 .3 1 from equation (j5]l with 
fc = 1.4-1.5 and a 10% uncertainty in the ratio /^jbsAfmodei- 
The limits on radii for the solar hydrogen composition are 
comparable to but a little lower than current theoretical expec- 
tations based on dense matter calculations which have radii 
of 10-13 km for neutron star equations of state that reach 
a maximum mass > 2Mq (Hebeler et al. 2010; Gandolfi, 
Carlson, & Reddy 2011). The mass-radius relation found in 
Steiner et al. (2010), derived from a set of photospheric ra- 
dius expansion X-ray bursts and hydrogen atmosphere fits for 
transiently accreting neutron stars in quiescence, also lies at 
slightly larger radii than our /?oo limits for solar composition. 
It should be noted that Suleimanov et al. (2011a) call into 
question the results of Steiner et al. (2010) by suggesting that 
"short" PRE bursts should be excluded from analysis as they 
show smaller blackbody normalizations in the burst tail and 
also do not follow the theoretically expected spectral evolu- 
tion. The implication is that the mass-radius relation derived 
in Steiner et al. (2010) would shift to higher radii as a result 
of using the more reliable "long" PRE bursts, and thus farther 
away from our derived upper limits. 

A smaller hydrogen fraction at the photosphere in 
GS 1826-24 would increase the R^o limits and make them 
consistent with the theoretical calculations and the mass- 



9 



radius curve from Steiner et al. (2010). We can get an im- 
pression of what the upper limit on Roo would be for an at- 
mosphere with a reduced hydrogen fraction by first looking at 
the extreme case of the pure helium atmosphere, and its de- 
rived upper limit of 21.5 km for logj,)(g) = 14.3 (see Figure 
E\. Such an upper limit is consistent with theoretical calcu- 
lations and the results from Steiner et al. (2010). We can go 
one step farther and estimate the hydrogen fraction we would 
require to have such a consistency with previous results using 
Equation [T9] Assuming a surface gravity of log[Q(g) = 14.3, 
an upper mnit on Rao of ~'16 km or more would be required. 
Using values for F^^^ and A averaged across the solar H/He 
model fits with log[Q(g) = 14.3, we estimate that a hydrogen 
fraction of X ss 0.5 or less would be needed. We are assum- 
ing that such a spectral model would not differ too greatly in 
shape from the solar H/He models. Galloway et al. (2004) 
find that the theoretical variations in burst properties with per- 
sistent flux between ignition models with X=0.7 and X=0.5 
are largely indistinguishable. However, more burst lightcurve 
simulations would be necessary to establish whether agree- 
ment with observed lightcurves is still possible with a reduced 
accreted hydrogen fraction. 

Given that the upper limit on Rrx, depends on the color cor- 
rection as f* (via A in equation 19 1, even a small 5% increase 
in the value of would yield a ^ 22% increase in the upper 
limit on Roo- Furthermore, Suleimanov et al. (2011b) dis- 
cuss large color correction factors of fc = 1.6-1.8 (cf. also 
Suleimanov & Poutanen 2006) possibly arising from a spread- 
ing layer associated with accretion onto the neutron star equa- 
tor. Perhaps in GS 1826-24 something similar is happening, 
although the increase in color correction required is not as 
large. It is worth noting that changing the visible area, for ex- 
ample by blocking one hemisphere of the neutron star with the 
accretion disk during the burst, does not change the inferred 
limits on radius because the limit on Roo is independent of 
the anisotropy factor (isotropic emission from only half the 
area is equivalent to setting = 2). 

There are several points to keep in mind when looking at 
our derived constraints on M and R. First, the constraints 
are only partly self-consistent in the sense that the lightcurve 
model used to fit the data does not have the same gravity as 
the derived M and R. Heger et al. (2007) (and Woosley et 
al. 2004) used a specific choice of gravity in their X-ray burst 
simulations. As we argue in §4, the lightcurve probably does 
not depend too sensitively on gravity, but additional simula- 
tions are needed to check this, and to calculate the uncertainty 
in the predicted model flux which enters in equation (|5]l re- 
lating fe and l+z- On the other hand, for the comparison to 
spectral models, in Figure |5] the upper limits on Roo, repre- 
sented by the solid coloredcurves, are placed in such a way 
as to coincide with the appropriate curve of constant surface 
gravity, consistent with the atmosphere spectral models used 
to derive those upper limits. 

A second issue is that the upper limit on Roo from the spec- 
tral models is based on fitting the initial part of the cooling tail 
only. We found that at first the slope of /T"'/^ with flux agrees 
well with the theoretical models of . In the latter part of the 
burst, however, at lower fluxes, the agreement breaks down. 
For F/FEdd ^ 0.2-0.3, /r"'/"* increases with decreasing flux, 
but more rapidly than expected based on the predicted val- 
ues, particularly those given by the solar metallicity models. 
Some other explanation is required for the rapid increase in fc 
and corresponding decrease in blackbody normalization in the 



tail of the burst. In 't Zand et al. (2009) suggest that this de- 
crease could be due to incorrect subtraction of the persistent 
emission, in particular, the subtraction of a thermal compo- 
nent that comes from the neutron star surface during accre- 
tion which is no longer present during the burst. Van Paradijs 
& Lewin (1986) pointed out that this effect should become 
important during the tail of the burst, when the burst flux be- 
comes comparable to that of the accretion. Looking at Figure 
[3] the observations and the low metallicity (solar metallicity) 
models begin to deviate at fluxes below ^ 5 x 10"^ erg cm"^ 
s"'(> 10"** erg cm"^ s~') compared to the persistent flux of 
2.1 X 10"^ erg cm"^ s"'. 

The disagreement between the observations and models be- 
gins sooner following the burst peak for solar metallicity than 
for low metallicity models because the former have a depres- 
sion in fc at low fluxes F /pEdd 0.3 (see Figure [3]), arising 
from absorption edges in partially-ionized Fe (Suleimanov et 
al. 2011b). There is no sign of such a dip in the observa- 
tions of GS 1826-24. This suggests a low metallicity in the 
photosphere, contrasting with the conclusions of Galloway et 
al. (2004) and Heger et al. (2007) who argued that the metal- 
licity was solar, based on burst lightcurves and energetics. A 
way to reconcile these disparate results is to consider the pos- 
sibility that Fe, whose presence has a significant influence on 
the spectrum but not on the burst energetics, may be absent 
from the atmosphere during bursts. If accretion halts during 
the burst, then Fe will rapidly sink through the atmosphere 
(Bildsten, Chang, Paerels 2003). On the other hand, since 
bursts from GS 1826-24 are all sub-Eddington, accretion may 
continue during the burst, resupplying Fe to the photosphere. 
Furthermore, while disk accretion only deposits mass near the 
equator, the accreted mass spreads faster latitudinally (< 0. 1 s; 
Inogamov & Sunyaev 1999; Piro & Bildsten 2007) than the 
timescale for Fe to sink through the atmosphere of ^ 1 s. Pro- 
ton spallation could also destroy a substantial amount of the 
accreted Fe (Bildsten, Chang, Paerels 2003). It should be 
noted that X-ray bursts produce a wide range of elements in 
their ashes which could significantly alter the spectrum. How- 
ever, mixing of burned material to the photosphere is thought 
not to occur due to the substantial entropy barrier (Joss 1977; 
Weinberg, Bildsten, & Schatz 2006). 

Another important unresolved issue is the changing spectral 
normalization K with accretion rate. The blackbody normal- 
ization K is 20% smaller for the 4.07 hr recurrence time 
bursts than the 5.74 hr recurrence time bursts. We cannot ex- 
plain this difference by changing the composition at the pho- 
tosphere and therefore changing fc (see discussion in §5.3). 
Also, it seems unlikely that a major change in composition 
would occur with only a w 50% change in accretion rate and 
smaller change in burst energy and lightcurves. As mentioned 
previously in §6, Suleimanov et al. (201 lb) discuss large color 
correction factors associated with accretion onto the neutron 
star equator, which they suggest accounts for the variations 
in measured K for the different spectral states of 4U 1724- 
307. Whether the 50% increase of accretion rate seen in GS 
1826-24 could result in the amount of hardening of the spec- 
trum observed needs to be investigated. If disk accretion onto 
the star significantly hardens the burst spectrum, it consider- 
ably complicates inference of the mass and radius from burst 
observations, and means that the range of fc included when 
calculating en^ors on mass radius determinations should allow 
for a larger range of values than given by spectral models. 

We have found that, even though they do not reach Ed- 
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dington luminosity, the bursts from GS 1826-24 show enough 
dynamic range in flux as they cool to significantly constrain 
^Edd by comparing with spectral models. A promising source 
to look at further is KS 1731-254 which shows both mixed 
H/He bursts with similar spectral evolution to GS 1826-24 
(Galloway & Lampe 2011) and photospheric radius expan- 
sion bursts. Analysis of these different bursts, which occur at 
different persistent fluxes and involve different fuel composi- 
tions (based on their energetics, peak luminosities and dura- 
tions), would give a stringent test of the spectral models and 
help to constrain any additional spectral components. 
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